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Abstract 

We consider within QCD collinear factorization the inclusive process p + p — >• 
hi + h2 + X, where the pair of identified hadrons, /ii,/i2) having large transverse 
momenta is produced in high-energy proton-proton collisions. In particular, we con- 
centrate on the kinematics where the two identified hadrons in the final state are 
separated by a large interval of rapidity Ay. In this case the (calculable) hard part of 
the reaction receives large higher order corrections ~ a" Ay". We provide a theoreti- 
cal input for the resummation of such contributions with next-to-leading logarithmic 
accuracy (NLA) in the BFKL approach. Specifically, we calculate in NLA the vertex 
(impact-factor) for the inclusive production of the identified hadron. This process has 
much in common with the widely discussed Mueller-Navelet jets production and can 
be also used to access the BFKL dynamics at proton colliders. Another application of 
the obtained identified-hadron vertex could be the NLA BFKL description of inclusive 
forward hadron production in DIS. 
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1 Introduction 



The BFKL approach [Ij is the most suitable framework for the theoretical description of the 
high-energy limit of hard or semi-hard processes, i.e. processes where a hard scale exists 
that justifies the application of perturbative QCD. 

At high squared center of mass energy s and for squared momentum transfered t not 
growing with s, large logarithms of the energy compensate the small coupling and must 
be resummed at all orders of the perturbative series. The BFKL approach provides a sys- 
tematic way to perform the resummation in the leading logarithmic approximation (LLA), 
which means resummation of all terms (asln(s))", and in the next-to-leading logarithmic 
approximation (NLA), which means resummation of all terms a<((as ln(s))". 

In the BFKL approach, both in the LLA and in the NLA, the high-energy scattering am- 
plitudes are expressed by a suitable factorization of a process- independent part, the Green's 
function for the interaction of two Reggeized gluons, and process-dependent terms, the so- 
called impact factors of the colliding particles (see, for instance, [2]). 

The Green's function is determined through the BFKL equation, which is an iterative 
integral equation, whose kernel is known at the next-to-leading order (NLO) both for forward 
scattering (i.e. for t = and color singlet in the t-channel) [U H] and for any fixed (not 
growing with energy) momentum transfer t and any possible two-gluon color state in the 
t-channel |5]. 

The impact factors of the colliding particle are a necessary ingredient for the complete 
description of a process in the BFKL approach and, therefore, to get a contact with phe- 
nomenology. The only impact factors calculated so far with NLO accuracy are those for 
colliding quark and gluons O [TJ |8], for forward jet production [9l [TOl [11], for the 7* — )■ 7* 
transition ^12j and for the 7* to light vector meson transition [13]. Recently, the NLO 
7* — )■ 7* impact factor has been calculated in the coordinate representation [T3|. 

The impact factors for colliding partons [6|, [7[ [8] played a role in the proof of fulfillment 
of bootstrap conditions for the gluon Reggeization in the NLA [2j and are at the basis of the 
calculation of the impact factors for forward jets. 

The impact factor for the 7* — )■ 7* transition is certainly the most important one from the 
phenomenological point of view, since it enters the prediction of the 7*7* total cross section, 
which is believed to be the golden channel for the manifestation of the BFKL dynamics |15j . 
Its calculation turned out to be very complicated and took several year to be completed |12j . 
A relatively simpler calculation led to the leading-twist NLO impact factor for the transition 
from a virtual photon 7* to a light neutral vector meson V = p'',a;,0, which was used to 
study the 7*7* — )■ VV reaction [16]. However, both for the 7*7* total cross section and for 
the production of two light vector mesons in 7*7* collisions, accurate experimental data will 
possibly come only from future high-energy e^e^ and 67 colliders. 

On the other side, the availability of the NLO forward jet impact factors [9j allowed 
the NLA calculation [T7] of the cross-section for the production of Mueller-Navelet jets [TS] . 
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which can be tested at hadron coUiders, such as Tevatron and LHC. 

In the present work we calculate in the NLO a new impact factor, namely the one for the 
forward production of an identified hadron, which can be used to determine the cross section 
for the semi-inclusive process p + p ^ hi + h2 + X , in the kinematics where the identified 
hadrons in the final state, hi, /i2, have large transverse momenta and are separated by a 
large interval of rapidity Ay. This process has much in common with the widely discussed 
Mueller-Navelet jets production and can be also used to access the BFKL dynamics at proton 
colliders. Very similarly to the case of forward jet production, the calculation of this impact 
factor lies at the border between coUinear and BFKL factorization: the proton in the initial 
state emits a parton which, in its turn, produces the hadron in the final state. This calls for 
the collinear factorization of the hard part of the process with parton distribution functions 
(PDFs) in the initial-state proton and with the relevant fragmentation functions (FFs) for 
the production of the final-state hadrons, both functions obeying the standard DGLAP 
evolution [T9] . 

The paper is organized as follows. In the next Section we will present the factorization 
structure of the cross section, recall the definition of BFKL impact factor and discuss the 
treatment of the divergences arising in the calculation; in Section 3 we give the derivation of 
LO impact factor; Section 4 is devoted to the calculation of the NLO impact factor; finally, 
in Section 5 we summarize our results. 



2 General framework 

We consider the process 

p{pi) + p{p2) ^ hi{ki) + h2{k2) + X (1) 

in the kinematical region where both identified hadrons have large transverse momenta^, 
fci ~ fc| 3> Aqqj-). This provides the hard scale, ~ fcf 25 which makes perturbative QCD 
methods applicable. Moreover, the energy of the proton collision is assumed to be much 
bigger than the hard scale, s = 2pi ■ p2 ^ k\2- 

We consider the leading behavior in the l/Q-expansion (leading twist approximation). 
With this accuracy one can neglect the masses of initial protons and identified hadrons. The 
state of the identified hadrons can be described completely by their (pseudo) rapidities 
and transverse momenta ki 2- 

In QCD collinear factorization the cross section of the process reads 

1 1 

—= / I dXidX2fi{Xi, jJ,)fj{x2, fi) — " {^1^2S, ^ ^2"^ 

dyidy2(Pki(Pk2 i^j=q,g-l dyidy2d'^kid'^k2 



^See Eq. ([3]) below for the definitfon of the transverse part of a 4- vector. 

^For massless particle the rapidity coincides with pseudorapidity, y — f], the latter being related to the 



particle polar scattering angle by 77 = — In tan 



2 • 
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where the i,j indices specify parton types, i,j = q,q,g, fi{x,^) denotes the initial protons 
parton density functions (PDFs), the longitudinal fractions of the partons involved in the 
hard subprocess are xi^2, is the factorization scale and da{xiX2S, ^) is the partonic cross 
section for the production of identified hadrons. The latter is expressed in terms of parton 
fragmentation functions (FFs), to be specified later. 

It is convenient to define the Sudakov decomposition for the identified-hadron momentum, 

P 

kh = othPi H -P2 + kh± , kl^ = -kl , (3) 

OihS 

where the longitudinal fraction an is related to the hadron rapidity as y = | In dy = ^ 

h ^ 

in the center of mass system. 

Another important longitudinal fraction is related with the coUinear fragmentation of the 
parton i into a hadron h. The probability for such fragmentation is generically expressed as 
the convolution of a parton fragmentation function and a coefficient function , 

1 

dai = C^{z)dz da^ = dan j ^D,^ /i) Cf (z, /x) . (4) 

Here is the cross section (calculable in QCD perturbation theory) for the production of a 
parton with momentum fraction z; the non-perturbative, large-distance part of the transition 
to a hadron with momentum fraction ah-, is described in terms of the fragmentation function 
D^] fi stands here for the factorization scale. 

Thus, the cross section da{xiX2S, /x) in Eq. ([2]) can be further presented as the convolution 
of two fragmentation functions with the cross section of the partonic subprocess 

i{xipi) +j{x2P2) m{ki/zi) +n{k2/z2) + X , (5) 

where X stands for inclusive production of additional, "unidentified" partons. This partonic 
cross section is free from infrared divergences, but it contains the coUinear singularities 
which are absorbed into the definition of PDFs and FFs, thus leading to a finite result for 
our process of interest ([1]) in coUinear factorization. 

In the case of large interval of rapidity between the two identified hadrons, the energy of 
the partonic subprocess is much larger than hadron transverse momenta, X1X2S ^ k"^. In this 
region the perturbative partonic cross section receives at higher orders large contributions 
~ a" In" related with large energy logarithms. It is the aim of this paper to elaborate the 
resummation of such enhanced contributions with NLA accuracy using the BFKL approach. 

Let us remind some generalities of the BFKL method. Due to the optical theorem, 
the cross section is related to the imaginary part of the forward proton-proton scattering 
amplitude, 

a = — . (6) 
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In the BFKL approach the kinematic hmit s ^ of the forward amphtude may be presented 
in D dimensions as follows: 

<5+joo 
<5— joo 

(7) 

where the Green's function obeys the BFKL equation 

c^G.(gi,g2) = 5''^'(gi-g2) + j d''~^qK{q,,q)GMqi) ■ (8) 

The energy scale parameter sq is arbitrary, the amplitude, indeed, does not depend on its 
choice within NLA accuracy due to the properties of NLO impact factors $1^2 to be discussed 
below. 

What remains to be calculated are the NLO impact factors $1 and $2 which describe the 
inclusive production of the identified hadrons hi and /12, with fixed transverse momenta ki, 
k2 and rapidities yi, y2, in the fragmentation regions of the colliding protons with momenta 
Pi and P2, respectively. For definiteness, we will consider the case when the identified hadron 
belongs to the fragmentation region of the proton with momentum pi, i.e. the hadron is 
produced in the collision of the proton with momentum pi off a Reggeon with incoming 
(transverse) momentum q. 

Technically, this is done using as starting point the definition of inclusive parton impact 
factor, given in Refs. [6l|7], for the cases of incoming quark(antiquark) and gluon, respectively. 
This definition involves the integration over all possible intermediate states appearing in the 
parton- Reggeon collision, see Fig. 1. Up to the next-to- leading order, this means that we 
can have one or two partons in the intermediate state. Here we review the important steps 
and give the formulae for the LO parton impact factors. 

Note that both the kernel of the equation for the BFKL Green's function and the parton 
impact factors can be expressed in terms of the gluon Regge trajectory, 

j{t) = 1 + u{t) , (9) 

and the effective vertices for the Reggeon-parton interaction. 

To be more specific, we will give below the formulae for the case of forward quark impact 
factor considered in D = 4 + 2e dimensions of dimensional regularization. We start with the 
LO, where the quark impact factors are given by 




where q is the Reggeon transverse momentum, and Tag denotes the Reggeon-quark vertices 
in the LO or Born approximation. The sum {a} is over all intermediate states a which 
contribute to the q ^ q transition. The phase space element dpa of the state a, consisting 
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of particles with momenta ^n■, is {pq is the initial quark momentum) 

while the remaining integration in flTOj) is over the squared invariant mass of the state a, 

Ml = iPq + qf . 

In the LO the only intermediate state which contributes is a one-quark state, {a} = q. 
The integration in Eq. ffTOj) with the known Reggeon-quark vertices T^qq is trivial and the 
quark impact factor reads 

= 9'^^^^ . (12) 
where g is the QCD coupling, as = g'^/lAir), = 3 is the number of QCD colors. 



In the NLO the expression ( fTOl) for the quark impact factor has to be changed in two 
ways. First, one has to take into account the radiative corrections to the vertices, 

rio) ^ F - r(°) + r(^) 

^ qq ^ ^ QQ ^ qq ^ ^ qq ■ 

Second, in the sum over {a} in (fTOj) . we have to include more complicated states which 
appear in the next order of perturbative theory. For the quark impact factor this is a state 
with an additional gluon, a = qg. However, the integral over becomes divergent when 
an extra gluon appears in the final state. The divergence arises because the gluon may be 
emitted not only in the fragmentation region of initial quark, but also in the central rapidity 
region. The contribution of the central region must be subtracted from the impact factor, 
since it is to be assigned in the BFKL approach to the Green's function. Therefore the result 
for the forward quark impact factor reads 

t2\ 



%{q, ^O) = 4 E / ) [raq{q)r ^(.A - M^] 

^ {a} 

_l/.-.|!*<o),J«o,,g,,-).,(^). ,13) 

The second term in the r.h.s. of Eq. f lT^ is the subtraction of the gluon emission in the 
central rapidity region. Note that, after this subtraction, the intermediate parameter s\ in 
the r.h.s. of Eq. (fT3|) should be sent to infinity. The dependence on sa vanishes because of 
the cancellation between the first and second terms, xi^^ is the part of LO BFKL kernel 
related to real gluon production. 

The factor in Eq. f[T3]l which involves the Regge trajectory arises from the change of energy 
scale (g^ Sq) in the vertices F. The trajectory function uj{t) can be taken here in the 
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xpi 




xpi 



Figure 1: Diagrammatic representation of the forward parton impact factor. 



xpi 




(Q/iPi, h 




Figure 2: Diagrammatic representation of the vertex for the identified hadron production 
for the case of incoming quark (left) or gluon (right). 



one-loop approximation {t = —q'^), 

N r d^-^k 



Uj{t) 



gH 



{2< 



D-l 



.^r(i-e)r-(e) 



(15) 



In the Eqs. f lTOj) and f|T3!) we suppress for shortness the color indices (for the explicit form 
of the vertices see [6j). The gluon impact factor ^g{q) is defined similarly. In the gluon case 
only the single-gluon intermediate state contributes in the LO, a = g, which results in 

^fl?-) =§;*r(^"')' (16) 

here Ca = N and Cp = (N'^ — 1)/{2N). In the NLO additional two-gluon, a = gg, and 
quark- ant iquark, a = qq, intermediate states have to be taken into account in the calculation 
of the gluon impact factor. 

The identified hadron production vertex we want to calculate is simply related with 
inclusive parton impact factors. In order to allow the inclusive production of a given hadron, 
one of these integrations in the definition of parton impact factors is "opened" (see Fig. 2). 
This means, in practice, that the integration over the momentum of one of the intermediate- 
state partons is replaced by the convolution with a suitable fragmentation function. We 
illustrate the procedure starting from the LO impact factor, contributing to the BFKL 
amplitude in the LLA, then we move on to the NLO impact factor, relevant for the BFKL 
resummation in the NLA. 
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3 Impact Factor in the LO 



The inclusive LO impact factors of quark (antiquark) and gluon are given in Eqs. (fT2!) 
and fll6p . It is important that these expressions are vahd also in non-integer dimensions. 

The inclusive LO impact factor of proton may be thought of as the convolution of quark 
and gluon impact factors with the corresponding proton PDFs, 



\ a=q,g / 

In order to establish the proper normalization for the impact factor in the case of frag- 
mentation with an identified hadron, we insert into the inclusive impact factor, Eq. ( IT7|) . a 
delta function and the integration over the "parent parton" transverse momentum, d'^k, and 
use Eq. (j4]) to get: 



^2 



a=q,q 



an expression for the LO impact factor schematically presented in Fig. |3l 

What is left is to express the "parent parton" variables by those of the identified hadron, 
k = {x/ah)kh (note that d^kh/kf^ = d'^k/k^). 

Although the results for the fragmentation impact factors should be expressed in terms 
of quantities describing the identified hadron, ah,kh, it is convenient during the calculation 
to operate with the kinematical variables of the parent parton, a,k (in LO, a = x). The 
transition to the hadron variables may be easily done at the end of the calculation. 

In what follows we will calculate the projection of the impact factor on the eigenfunctions 
of LO BFKL kernel, i.e. the impact factor in the so called [u, ra) -represent at ion, 

^u,n) = [ d'q^^ (g2) — ^e-<^ . (19) 
J q 7rv2 

Here (f) is the azimuthal angle of the vector q counted from some fixed direction in the 
transverse space. 



4 The NLO calculation 



We will work m D = A + 2e dimensions and calculate the NLO impact factor directly in the 
(z/, n)-representation flTIJl) . working out separately virtual corrections and real emissions. To 
this purpose we introduce the "continuation" of the LO BFKL eigenfunctions to non-integer 
dimensions, 

(g^)^e-^^(g^)^-5(g-r)" , (20) 
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Figure 3: Diagrammatic representation of the LO vertex for the identified hadron production 
for the case of incoming quark (left) and gluon (right). 

where 7 = zz^ — | and / ^ = 0. It is assumed that the vector I lies only in the first two of the 

2 + 2e transverse space dimensions, i.e. I = ei + i 62, with efg = 1; el " 62 = 0. In the limit 
e — 7- the r.h.s. of Eq. fj20|) reduces to the LO BFKL eigenfunction. This technique was 
used recently in Ref . [20j . An even more general method, based on an expansion in traceless 
products, was uses earlier in Ref. [21] for the calculation of NLO BFKL kernel eigenvalues. 



Thus, for the case of non-integer dimension the LO result for the impact factor reads 



^TTT = C dah 



which gives in the (z/, n)-representation the result 



(21) 



n) 



/^(|/,w^;(^).E/.w-;(f))(* 



k-l 



Oh 



a=q,q 



(22) 



Collinear singularities which appear in the NLO calculation are removed by the renor- 
malization of PDFs and FFs. The relations between the bare and renormalized quantities 
are 

f,{x) = f,ix, /.^) - t (I + In ^) / f [P,,(^)/,(f , /ii.) + P..(^)/.(f , I^f)] , 

X 

f,{x) = /,(x,/x^) - ^ (1 + Inf ) / f [P,,(^)/,(|,/xp) + P,,(^)/,(f ,/i^)] , (23) 



where | = 7 + 7_b — ln(47r) ^ ^l^^w , and the DGLAP kernels are given by 



C, 



e(47r) 



Tn [z' + (1 - zf] 
2Ca ^ 



1 + '^' 3^,^ 



(24) 
(25) 

(26) 

(27) 



with Tr = 1/2. Here and below we always adopt the MS scheme. Similarly, for FFs we have 
D':^{x) = D^{x,fiF) - t (i +ln^) /f K(f,Mn.(^) +^.'(f,/"F)n,W] , 

X 

D^{x) = D^^ix, lip) - ^ (4 + In f ) / f [D^^{% fiF)P,,{z) + I),^(f , liF)P,,{z))] • (28) 



Here /i is an arbitrary scale parameter introduced by the dimensional regularization, which 
cancels out in the results for physical quantities. Note the difference in the non-diagonal 
terms {Pgg Pgq) between the formulas for PDFs' and FFs' renormalization, which is due 
to the fact that in the parton to hadron splitting the hadron is in the first case in the initial 
state, whereas in the second case it is in the final state. 



Owing to k = {x/ah)kh, one gets 



(29) 



k^]' M k,, ■ I 




Now we can calculate the collinear counterterms which appear due to the renormalization 
of bare PDFs and FFs. Inserting the expressions given in Eqs. ( 123|) . (128!) into the LO impact 
factor (l29l) . we obtain, after some suitable transformations of the integration variables in the 
terms coming from the PDFs renormalization, the following collinear countertermEl 



7tV2P d^''{i^,n)\coii. c.t. 



C 



dahd^+^'k 



27C 



X 



z 



k-l 



Oh 



:i + z 



'27- 



P„(r.) (^) + (§i + r.-A P„(z) J: 

a=q,q ^ ' <i=<?i9 



(30) 



The other counterterm is related with QCD charge renormalization; with nf active quark 
flavors we have 



as = asifiR) 



l + ^/3of'+ln4 



4tt 



_ IICa 2nj 

; PO — n ^ ' 



(31) 



■^In principle, one can consider different values of the factorization scale for the evolution of PDFs and 
FFs. Here for simplicity we take these scales equal. 
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and is given by 



1 

X 
ah 



/ ? (§;«^)-» (f) - E (^)) (-0" • (-) 



a=q,q 



To simplify formulae, from now on we put the arbitrary scale of dimensional regularization 
equal to the unity, /i = 1. 

In what follows we will present intermediate results always for ^v^^^ 1'^'?'*^ , which we 
denote for shortness as 

Note that we always denote by k the transverse momentum of the parton which fragments 
to the identified hadron. Moreover, with no argument from now on is to be understood 
as asifiR). 

We will consider separately the subprocesses initiated by a quark and a gluon PDF, and 
denote 

/ = /, + /,. (34) 
We start with the case of incoming quark. 



4.1 Incoming quark 

We distinguish virtual corrections and real emission contributions, 

I, = + ■ (35) 



Virtual corrections (see Fig. H]) are the same as in the case of the inclusive quark impact 
factor, therefore we have 



a=q,n 



xlCpi- — + A - rif- ^ + 4^ + ^(1 - e) - 2^(e) + ^(1) 

^ l + 2e J ^(l + 2e)(3 + 2e) ^\ ^\ ) ^\ ) 



^4(l + 2e)(3 + 2e) e(l + 2e) 4(1 + 2e) 2^ ' ^^^^ 
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+ 



Figure 4: Diagrammatic representation of the NLO vertex for the identified hadron produc- 
tion for the case of incoming quark: virtual corrections (the larger black blobs indicate the 
1-loop contribution). 

We expand fl36ll in e and present the result as a sum of the singular and the finite parts. 
The singular contribution reads 

271 {AttY eT{l + 2e)J x ^ ' " \ x ) \ ) V / 

ah " '^''^ 

2 ^/ , ^ A -^o , 11 



x^CH;-3l-f + C.ll„^ + -|^ (37) 



whereas for the regular part we obtain 



1 



a=q,q 



5nf _ /85 . TT^ 



^<^C, + ^-Ca{- + -]} . (38) 



Note that {l^)^ + differs from by terms which are 0{e). 



4.1.1 Quark-gluon intermediate state 

The starting point here is the quark-gluon intermediate state contribution to the inclusive 
quark impact factor, 

- ^.^V'gg^f ^tjtg,. {C^^i ^ 0^,. - ,A . ,) ] . (30, 

— * — * 

where /5i and are the relative longitudinal momenta (/3i + P2 = 1) and ki and ^2 are the 
transverse momenta (fci + k2 = q) oi the produced gluon and quark, respectively. 

We need to consider separately the fragmentation of the quark and of the gluon into the 
final hadron, therefore we denote the corresponding contributions as I^^ and 7^^, 

Q m ^ q,9 ^ 
We start with the gluon fragmentation case, I^ g (see Fig. |5]). 
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Figure 5: Diagrammatic representation of the NLO vertex for the identified hadron produc- 
tion for the case of incoming quark: real corrections from quark-gluon intermediate state, 
case of gluon fragmentation. 

a) gluon fragmentation 

The "parent parton" variables are k = ki, ( = f3i {(32 = C = 1 — C; ^2 = Q — k), therefore 
we have 



1 1 



q,9 



a. 



2n{ATcy J TT 



[q ) M g ■ / ' 



h 0!h 



X J C \xC 



^ - k- q 



kfiq- 



(41) 



It is worth stressing the difference between the previous calculations of NLO inclusive par- 
ton impact factors and the present case of fragmentation to a hadron with fixed momentum. 
In parton impact factor case, one keeps fixed the Reggeon transverse momentum q and inte- 
grates over the allowed phase space of the produced partons, i.e. the integration is of the form 
/ 2C(i-c) • • -'■^ fragmentation case, instead, we keep fixed the momentum of the 

parent parton /c, and allow the Reggeon momentum q to vary. Indeed, the expression ( HTll 
contains the explicit integration over the momentum q with the LO BFKL eigenfunctions, 
which is needed in order to obtain the impact factor in the (z/, n)-representation. 

The g-integration in ( 14T|) generates 1/e poles due to the integrand singularities aX q k/C, 
for the contribution proportional to Cp and at g — )■ for the one proportional to Ca- 
Accordingly we split the result of the g-integration into the sum of two terms: "singular" 
and "non-singular" parts. The non-singular part is defined as 
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a.. 



1 1 
dx f d( 



ah 



a=q,q 



c V c 



X 



— n ■ a 



1+e 



(a - n)2 ( a - i 



2\7- 



a ■ / 
n ■ I 



where n is a unit vector, = 1. Taking this expression for e = we have 



(42) 



2vr / X / C 
where we define the function 

h = A (^,7,0 



cv c 



h , (43) 



TT 



a — n r a — 



(44) 



with the azimuth of the vector a (0) counted from the direction of the unit vector n. 
For the singular contribution we obtain 



1 1 



ZTT (47r)^ e r(l + 2e) y a; J C ' ^ 



ah 



a=q,q 



k-l 



X 



l + C' + eC^ 



^ r(i + 26)r(f - 7 - 6)r(f + 1 + 7 + 6) /c 

^r(l + e)r(l-e)r(f-7)r(f + l + 7 + 2e)'' ""V 



2e 



Expanding it in e we get 



rvi, a=q,q 



k-l 



1 + C 



+ e 



2 r 



c 



CFr'"(x(^,7)-21nC) + 2C^ln 



+ C(C^Fr'" + C^)) 



where 



x(n,7) = 2V'(l)-V'(^-7) -V^(^ + l + 7 
is the eigenvalue of the LO BFKL kernel, up to the factor Nas/ir. 



(45) 



(46) 



Note that the divergent part of is canceled by the corresponding term of the coUinear 
counterterm ( 15U]) . 
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Figure 6: Diagrammatic representation of the NLO vertex for the identified hadron produc- 
tion for the case of incoming quark: real corrections from quark-gluon intermediate state, 
case of quark fragmentation. 



b) quark fragmentation 

Now the "parent parton" variables are k = C = {Pi = C, ki = q — k) (see Fig. EJ. 
The corresponding contribution reads 



a. 



27r(47r) 



TT 



q-l 



1 1 
dx f dQ 



X 

ah ^ 



C - V XL 

a=q,g ^ 



Oh 



X 



(1-0 



k'- 



(g- A;)2 (q 



Ca- 



k-q"-±< + f 



'-kyiq 



(47) 



We will consider separately the contributions proportional to and Ca- 
bi) quark fragmentation: CiT'-term 

Note that the integrand of the Ci?-term is not singular at — )■ 1. We use the decompo- 
sition 

/ \ 

1 1 

+ 



- ky + (q- I 



{q- kf 



q- f 



V J 



{q- [q 

in order to separate the regular and singular contributions. The regular part is given by 
1 1 



a. 



27r(47r)^ 



Oh 



ay, a=q,q 



c 



X Cf 



d^+^'d 



1+e 



(a-n)2 + (a- i 



fi-l 



(a — n)" 



1 id 
- + 



n-l 



[a — 



n\2 



(48) 
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Therefore for e = we have 
1 1 



a=q,q 



(49) 



where we define the function 

T T ( ^\ ■ 

h = /2(n,7,C) 



^ ia — nY + 



(a — fiy 



+ 



(50) 



The singular part is proportional to the integral 



/ 



{q- kf 



T{l-e)T\l + e) (^^V^^- 



therefore, for the singular part of the CF-term we have 

1 1 



2€-2 



r(l-e) r'(l + e) [ dx f dC 
2^ e(47r)^ r(l + 2e) 



ah o,=q,q 



X Ci 



1 + C + /C 



1-C 



(9^1 ^c-) 



The next step is to introduce the plus-prescription, which is defined as 



1 

/ 



no 

(1-0- 



1 

/ 



^(C) - i^(i) 
(1-C) 



a 

-I 



dC 



(1-C) ' 



(51) 



(52) 



for any function F{Q, regular at ^ = 1. Note that 

(1 - 0- ^ (1 - 01-' + ^^(1 - C) ^ ^^(1 - C) + + 2e (Mi-^ 



Using this result, one can write 



1-C VC 



2e 



I - 3 j S{1 - C) + [-^^^ + ^^(1 - ) (1 + C-'^) + 0{e) 
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= _ 3^ 5(1 _ ^) + p^^(^) (1 + ^-27) + o{e) . 

Taking this into account and expanding in e the singular part of the Ci?-term, one gets the 
following result for the divergent contribution: 



Cf _ as T[l - e] r2(l + e) f dx f dC 



X 



+ eCAl + C'-') ( C + 2(1 + C^) (^^) - 2(1 + C^)- 



1-C 



1-C 



k-l 



(53) 



Note that the first term of the divergent contribution given in ( l53l) cancels in the sum with 
the corresponding term of the virtual correction ( l37I) . whereas the second one is canceled by 
the corresponding term of the collinear counterterm (l30|l . 

quark fragmentation: C^-term 

The C^-contribution needs a special treatment due to the behavior of ( H7|) in the region 



We use the following decomposition 

Ca 



C) 



Ca 



2(1-0 



1 1 



2 (q-k)2{l-C) 



The second term in the r.h.s. is regular for C — )■ 1 and can be treated similarly to what we 
did above in the case of the CiT'-contribution. The first term is singular and the integration 
over C, has to be restricted (according to definition of NLO impact factor) by the requirement 



and assuming sa — oo. Therefore the ^-integral has the form 

i-Co 



dC 



no 

1-C 



for 



Using the plus-prescription ( 1521) one can write 

i-Co 1 



FjC) 

1-C 



d( .f^^), +F(l)lni for Co ^ 



(54) 
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for any function not singular in the limit C ~^ 1, and 



:i-c) 



Ca 



2 {q-kf {q-kf 



+ 



2 (^-_fc)2(l-0+ 2(1-0 



k\2 



(q-ky 



CJ (g-'_A;)2(g-_^)2 



We remind that the definition of NLO impact factor requires the subtraction of the con- 
tribution coming from the gluon emission in the central rapidity region, given by the last 
term in Eq. f ll3p . which we call below "BFKL subtraction term". After this subtraction the 
parameter s\ should be sent to infinity, sa — )■ oo. Our simple treatment of the invariant mass 
constraint, Mqq < sa, anticipates this limit s\ — ?■ oo, therefore we neglect all contributions 
which are suppressed by powers of l/s\. Moreover, the first term in the r.h.s. of the above 
equation should be naturally combined with the BFKL subtraction term, giving finally 







q^-k-q^^ + ^. 
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{q-ky [q 



+ 



+ 



2 (g-._^)2(l_C)^ 2(1-0 



^^5(l-0^^1n^^ 
2 [q-ky {q-ky 



{q- ky Vc 



k\2 



ky{q- f )2 



After that, we are ready to perform the integration, which naturally introduces the sepa- 
ration into singular and non-singular contributions. The singular contribution reads 



a,r(i-6) r^(i + 6) 

271 e(47r)^ 1(1 + 2e) 
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X ■ 
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■So 
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(5(1-0 ln^ + ^ --7-e +^ U + 7 + 7T + 2e -V^(e)-V^(l) 



n 



(C 



c 



-2t 



-27-2e> 



Expanding this expression in e and using that 



and 



1-C 



^2e-l ^^-2e _^ ^-27-2e^ ^ (^^-2e _^ ^-2-y-2e^ 



2e 



+ 
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+ 0{e 



;i-c)- 



(55) 
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Figure 7: Diagrammatic representation of the NLO vertex for the identified hadron produc- 
tion for the case of incoming gluon: virtual corrections (the larger black blobs indicate the 
1-loop contribution). 



we get a divergent term, 



Ca C(s r(l — e) f dx f d( 



a=q,q 



x<iC^5(l-C)ln| 



(56) 



+ eCA 



£_o , 1 

k 

2^ / n , .-2,^ ( X(n,7) _ / ln(l-C) \ \ InC 

2(1-0+ V 1-c ;J 1-c 



+ (l + a (1 + C 



The divergent term given in fISBl) cancels in the sum with the corresponding term of 
the virtual correction fl57|) . The remaining singularity in f l?r|) vanishes after the charge 
renormalization fl?I]) . 

The regular contribution differs from ( HHj) only by one factor and reads 



\ 1,1) r 



Ca Cis f dx f d( 



2vr J X J C 



h 

xC 



k-l 



(1 + C^) 



(57) 



4.2 Incoming gluon 

We distinguish virtual corrections and real emission contributions, 

h = Iq + If ■ (58) 

Virtual corrections (see Fig. [7]) are the same as in the case of inclusive gluon impact 
factor, 

^ 271 {Any eT{l + 2e) J X ■'^^ ' ^ \ X J \ J V / Cp 

a-h 
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Figure 8: Diagrammatic representation of the NLO vertex for the identified hadron pro- 
duction for the case of incoming gluon: real corrections from quark- ant iquark intermediate 
state, case of quark fragmentation. 

- n - 2(1 + 2el(3; 2e) + " " + + ^^^^ 

, 6 \ / (l + e)(2 + e)-l-^ \) 

^(l + e)(l + 2e)(3 + 2e)y' ^""^ V (l + e)(l + 2e)(3 + 2e) j ■ ^ 

Expanding it in e we obtain the following results for the singular, 

(jv. _ «. r[i-6] i r^(i + 6) I u^jjH (j:2y^^-"^ (j^^iY Ca 

>s - 2vr {Any e 1(1 + 2e)J x ^'^""^^^ U J T J ^ J 

Sq 2 11\ rif 



and the finite parts. 



277 / X ' ^ V a; / V ; V J C 



x\CA{'^^-i]-^.nA . (61) 



18 2 J 9 

For the corrections due to real emissions, one has to consider quark-antiquark and two- 
gluon intermediate states, 

^ = ^f, + • (62) 



4.2.1 Quark-antiquark intermediate state 



The starting point here is the quark-antiquark intermediate state contribution to the inclusive 
gluon impact factor (Tr = 1/2), 



CAk^ki k^ki{k2/3i-ki/32y 



(63) 
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where /3i and /32 are the relative longitudinal momenta + /32 = 1) and ki and k2 are the 
transverse momenta {ki + k2 = q ) of the produced quark and antiquark, respectively. 

We consider quark fragmentation (see Fig. |8]). In this case, the "parent parton" variables 
are k = ki, ( = /3i {(32 = ( = ! — (, k2 = q — k). The case when antiquark fragments is essen- 
tially the same. Therefore, the sum over all possible quarks and antiquarks fragmentation 
has the form 



9,1 



a. 



27r(47r) 



1+e 



ah ^ ""'^''^ 



X ( 1 - 



2CC\ Cp 1 , C k-iq-k) 



1 + e J [CA{q-ky C(g-^)2(g-^)2 



(64) 



We can split this integral into the sum of singular and non-singular parts. The non- 
singular contribution of (Elj) reads 



27r(47r) 



7+e- 



TT 



1+e 



{ay-- ^ 

\ fi ■ I 



Expanding it in e we obtain 



-27 



a ■ n — 1 



(d-ny (a- f 



2 • 



(65) 



ah ^ a-?."? 



^*"*^*-r)"^jp„(c)/3. (66) 



where we define the function 

h = /3(ri,7,C) 



(Pa a ■ n — 1 



TT 



(67) 



For the singular contribution of flMj) we have 



27r {AttY er(l + 2e) 



ah 



^ a=q,q ^ ^ 



xTfi 1 



2CC 
1 + e 



Cf r(i + 26)r(f-7-6)r(f + 1 + 7 + 6) 

[Ca r(l + e)r(l - e)r(f - 7)r(f + l + 7 + 2e) 



2e A-2e-27 



• (68) 
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Figure 9: Diagrammatic representation of the NLO vertex for the identified hadron produc- 
tion for the case of incoming gluon: real corrections from the two-gluon intermediate state, 
case of gluon fragmentation. 



Expanding it in e, we get 

\ 9,q/s 



271 e(47r 
Cf 



" Ca 
Cf 



+ c 



-27 



(69) 



+e(2CCT. (| + C-)+P.(C) (|x(^>7) + 2C-ln0)} . 



Note that the divergent part of this expression is canceled by the corresponding term of 
the collinear counterterm fl30l) . 



4.2.2 Two-gluon intermediate state 



The starting point here is the gluon-gluon intermediate state contribution to the inclusive 
gluon impact factor, 



PI 



+ 



/3, 



(70) 



where j3i and are the relative longitudinal momenta + /32 = 1) and ki and ^2 are the 
transverse momenta {ki + k2 = q) oi the two produced gluons. 

We need to consider the case of a gluon which fragments, the case when the other gluon 
fragments being taken into account by a factor 2 (see Fig. |9]). Thus, we obtain the following 
integral: 
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{q- kf (q 



(71) 



The calculation goes along the same lines as in the Section [4.1.11 (case b2). First, we 
separate the ( 1 singularity, then we add the BFKL subtraction term. Using flMjl one 
obtains 
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Again, we can split the result into the sum of singular and non-singular parts. The 
non-singular contribution differs from ( HHj) only by a factor, it reads 
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(72) 



Expanding it in e we obtain 
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(73) 



For the singular contribution we obtain 
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Expanding this result in e, we get for the divergent contribution 
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Finally, the e expansion of the divergent part has the form 

1 1 
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+2(1 + C 



(76) 



The term in the divergent contribution (!76l) proportional to Pgg cancels with the correspond- 
ing term of the collinear counterterm (150]) . while the term in f l76p proportional to 5(1 — () 
cancels in the sum with the singular part of the virtual corrections (1601) . The uncanceled 
divergence in the virtual corrections f l60|) vanishes after the QCD charge renormalization (|3T]) . 

Summarizing, all the infrared and ultraviolet divergences arisen in the calculation have 
disappeared after taking into account PDFs and FFs renormalization and QCD charge renor- 
malization. Collecting all intermediate contributions we obtain the final result for the iden- 
tified hadron NLO impact factor, which reads 
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The results for the NLO coefficient functions read 
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c,,(x,c) = P,,(C) 
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For the /i 2 3 functions we obtain the following results: 
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Using the following property of the hypergeometric function, 

2Fi(l,a,a + l,0 =a(^(l)-^(0-lnO + (^(ClnO , 
one can easily see that 

/2 = O(ln0, /i = O(ln0, /3 = O(ln0, 
which implies that the integral over C, in f l77|) is convergent in the upper limit. 
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5 Summary 



In this paper we have calculated the NLO vertex (impact factor) for the forward production 
of an identified hadron from an incoming quark or gluon, emitted by a proton. This is a 
necessary ingredient for the calculation of the hard inclusive production of a pair of rapidity- 
separated identified hadrons in proton collisions ([1]). This process, similarly to the production 
of Mueller- Navelet jets, can be studied at the LHC hadron collider. 

Another natural application of the obtained identified hadron production vertex could 
be the NLA BFKL description of inclusive forward hadron production process in DIS, 

e{pi) + p{p2) ^ h{k) + X , (85) 

where in the low- a; event the hadron h{k) with high transverse momentum is detected in the 
fragmentation region of incoming proton p{p2)- Data for such reaction in the case of forward 
7r°-production were published by the HI collaboration at HERA \22\ . 

At the basis of our calculation of the hard part of the vertex was the definition of NLO 
BFKL parton impact factors; then the collinear factorization with the PDFs of the incoming 
partons and with the FF for the production of the identified hadron (in the MS scheme) was 
suitably considered. 

We have presented our result for the vertex in the so called {u, n)-representation, which 
is the most convenient one in view of the numerical determination of the cross section for 
the production of a pair of rapidity-separated identified hadrons along the same lines as in 
Ref. [ig. 

We have explicitly verified that soft and virtual infrared divergences cancel each other, 
whereas the infrared collinear ones are compensated by the PDFs' and FFs' renormalization 
counterterms, the remaining ultraviolet divergences being taken care of by the renormaliza- 
tion of the QCD coupling. 

In our approach the energy scale Sq is an arbitrary parameter, that need not be fixed 
at any definite scale. The dependence on sq will disappear in the next-to-leading logarith- 
mic approximation in any physical cross section in which the identified hadron production 
vertices are used. Indeed, our result for the NLO vertex, given by Eqs. ( !77|) -( !8T|) . contains 
contributions ~ In(so) and these terms are proportional to the LO quark and gluon vertices 
multiplied by the BFKL kernel eigenvalue x('^) ^)- This fact guarantees the independence of 
the identified hadrons ([1]) or single hadron ( 185|) cross section on sq within the next-to-leading 
logarithmic approximation. However, the dependence on this energy scale will survive in 
terms beyond this approximation and will provide a parameter to be optimized with the 
method adopted in Refs. |16] . 
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A Useful integrals 



We give here some useful integrals: 



A;'2 \k'^ + {k'-kyj 2y ^/2(^/_^)2 V J er(l + 2e) 
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In the integrals below, = is assumed: 
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